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{DJO, The cylindrically-symmetric static vacuum equations of Conformal Gravity are solved for 

the case of additional boost symmetry along the axis. We present the complete family of 
solutions which describe the exterior gravitational field of line sources in Conformal Gravity. 
' We also analyze the null geodesies in these spaces. 
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1 Introduction 



^ \ Conformal Gravity (CG) [1] was proposed as a possible alternative to Einstein gravity ("GR"^ 



which may supply the proper framework for a solution to some of the most annoying problems of 
theoretical physics like those of the cosmological constant, the dark matter and the dark energy. 

The gravitational field in CG is still minimally coupled to matter, but the dynamical basis 
is different: it is obtained by replacing the Einstein-Hilbert action with the Weyl action based 
on the Weyl (or conformal) tensor Ci^\^y defined as the totally traceless part of the Riemann 
tensor: 

1 R 
CKXfMu = RkXijlu — -^{gKfiRxiy — QkuRxij. + gxuRn^i — gxtiRKu) + —{gnfiOxu — gKvgx^i), (l-l) 



so the gravitational Lagrangian is 



/:, = -^c,am.c"^^''" (1.2) 



*Electronic addresses: verbin@openu.ac.il; Yves.Brihaye@umons.ac.be 



where o; is a dimensionless parameter. The gravitational field equations take the following form: 

a 

W^u = ^T^. (1.3) 
where T^;^ is the energy-momentum tensor and W^i, is the Bach tensor given by 

W"" = R^xC^^^" - 2V xC^^'' , (1.4) 
or in terms of the Riemann and Ricci components by: 
W^, = ^V^V.i? - VAV^i?^. + ^(^' + ^AV^i? - 3R^xR^^)9^. + 2R^^Ri,^,x - ^i^i?/..- (1-5) 

It was suggested (see [1] and references therein) that while CG agrees with Newtonian grav- 
ity in Solar System scales, it further produces a linearly growing potential that could explain 
galactic rotation curves without invoking dark matter. It was further argued that accelerating 
cosmological solutions of CG describe naturally the accelerated expansion of the universe thus 
removing the need for dark energy. 

On the other hand, CG has been criticized from several aspects both phenomenological and 
formal. Arguments in favor of the need of dark matter come from observations of the unusual 
object called "bullet cluster" [2, 3] whose dynamics seems very difficult to understand without 
assuming a weekly interacting dark component. 

More specifically, several authors claim that predictions in the weak field limit of CG disagree 
with solar system observations [4] , yield wrong light deflection [5] and that the exterior solutions 
cannot be matched to any source with a "reasonable" mass distribution [6] . Other authors find 
evidence for tachyons or ghosts [7] or raise the fact that only null geodesies are physically 
meaningful in this theory since the "standard" point particle Lagrangian is not conformally- 
invariant [8, 9]. 

Counter arguments to some of these objections were also published [10, 11], as well as possible 
ways [9] out of some of the difficulties and the matter is, to our view, still waiting for a consensus. 

It is therefore very much required to investigate further the predictions and consequences 
of CG in its purely tensorial formulation as well as in its scalar tensor extension as much as 
possible. 

In this work we concentrate on cylindrically-symmetric static vacuum solutions with the aim 
of clarifying further the properties of string-like solutions in CG [12]. Cosmic strings [13] are 
a typical outcome in any field theory which describes matter in the very early universe, thus 
serving very well the purpose of testing the implications of CG. The main result reported here 
is the full family of the static cylindrically-symmetric vacuum solutions of CG which represent 
the external gravitational field of localized line sources. 
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2 Cylindrically-Symmetric Equations 



The general static cylindrically symmetric line-element has the form: 



(2.1) 



In order to find solutions for this system, we have to fix the arbitrariness of the radial 
coordinate and the arbitrary rescaling of the metric due to the conformal symmetry. Wc will 
also limit our solutions to those exhibiting boost symmetry along the string direction (z), as 
for ordinary cosmic strings, i.e. K{r) = B{r). This leaves one independent metric component. 
So only one of the gravitational field equations (1.3) (with Ty = 0) has to be solved. We may 
therefore choose to solve the lower order rr equation, i.e. = 0. 

Several special solutions are already known in explicit form. First of all, all the thin string 
(line source) solutions of GR with either a vanishing or non- vanishing cosmological constant [14] 
which solve for r > 



satisfy also the CG vacuum equations Wfj^v = 0. This is obvious by direct substitution of (2.2) 
inEq. (1.5). 

One special member of the A < family is the AdS soliton [15] (see also [16, 17, 18]) which 
is a cylindrically-symmetric regular solution of the same equation (2.2) and is therefore distinct 
from AdS (anti-de Sitter) space. More recently, we have discovered two families of very simple 
exact solutions [12] during a mainly numerical study of cylindrical solutions of the Abelian Higgs 
model coupled to CG. However, we were unable to integrate the equations analytically at the 
time. 

Here we present a reduction of the vacuum equations for static cylindrically-symmetric so- 
lutions to a single first order non-linear equation which may be solved by a straightforward 

quadrature. 

We have therefore first to complete the gauge fixing. Since we would like it to be consistent 
with the symmetric vacuum solutions [i.e. (A)dS spaces] and with the AdS soliton, it should 
respect the asymptotic condition 



where k is a real parameter. Note however, that a constant Ricci scalar is not a "gauge invariant" 
concept in CG; it is only a matter of convenience which can be obtained by a proper gauge choice. 
The corresponding invariant condition is that the Weyl tensor will vanish asymptotically, that is 
spacetime becomes conformally flat asymptotically. These restrictions will simplify considerably 



(2.2) 



^(5^ for r -)• oo 



(2.3) 
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the very cumbersome expressions of the components of the Bach tensor and will enable a clear 
physical picture. 

The simplest gauge choice is B{r) = K{r) = M{r) = 1 with L{r) as a single metric compo- 
nent. This choice may be obtained directly from Eq. (2.1) by a suitable conformal transformation 
combined with a redefinition of the radial coordinate. In this gauge Rq = = soii cannot 
contain the well-known solutions with a cosmological constant of Eq. (2.2). The cylindrical 
version of the "Mannheim gauge" [1] M(r) = 1/B{r) ,L{r) = r with K{r) = B{r), turns out to 
be complicated and the simplifications from spherical symmetry are not observed. 

An alternative parametrization of the metric tensor which solves these difficulties is the 
following gauge choice 

M = 1 , L = — = B' , Kir) = B(r) (2.4) 
dr 

where we also resort to dimensionless coordinates. This metric is equivalent to that of the 
"Mannheim gauge" , but has the advantage that it leads to autonomous equations which we can 
solve by quadrature. 

Note however, that this gauge excludes the "flat" A = solutions. Conformally flat solutions 
are of course still allowed. 

In this gauge the components of Ricci tensor and scalar take the form 



B'\^ B" B'" B" 

R^^ = Rl = -[—\ -2— , r; = R'P = -_-2— , R 



B'V B" B'" 

b) 5^ 



(2.5) 



while the rest vanish. 

Already at this stage we can obtain easily the constant Ricci solutions of Eq. (2.2). It is 
enough to solve the (00) equation which is readily integrated to give 

\B'f + —B^--=0 (2.6) 
2^ ^ 24 5 ^ ^ 

where b is the integration constant. This is a trivial "mechanical" equation whose solutions are 
easily obtained by inspection. We will discuss these solutions within the more general framework 
in the next section. 

The non-zero components of the Weyl tensor are all proportional to a single quantity namely 
C%, = C\=C^\, = C'^%, = C/6 , C^\=C'\, = -C/3 (2.7) 



where 



^see the discussion about the metric g^iv = diag{l, —1, —H^, —1) at the end of this section 
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and its square is given simply by C^.XfxvC'^^^'^ = 4C^/3. 
the non-vanishing components of the Bach tensor are 



5(5) 

+ 



3S' 



35 



+ 



2B"'B" 2B"'B' 



SB'"^ 3B'^ ' 3BB' 
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(2.9) 
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(2.10) 



while the fourth one, W!p can be obtained immediately from the identity Wf^ = 



Using an exponential transformation B 



al3 



simplifies the expressions a little. VF^ and the 



"Weyl tensor quantity" C which will be used in the following, become: 

)2 



= - 



„2 ^ 2/3(3) ^"2 4^(3)^// 2^(4) /3" 
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(2.11) 



3/?'2 3/3'3 3/3' 3/3'2 ' 

An easy way to obtain is to start with the line clement ds^ = B^ \dt^ — dg^ — H'^dip^ — dz^\ 
whose Bach tensor is given by = W^/B'^ where is calculated from the metric g^j^n = 
diag{l, —1, —H"^, —1). The components Wj^ are relatively easy to obtain since ^^j, has only one 
non- vanishing independent Riemann component: Rg^g^ = —HdpH/dg^. We thus get 
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IF' 



A. 
dg 



(2.12) 



while W!p can be obtained from Wl^ 



0. A further coordinate transformation Bdg = dr 
together with H = d(3/dr yields after lengthy calculations the above components of the Bach 
tensor of Eqs. (2.9)-(2.10). As a check we have verified that these components satisfy a covariant 
conservation law, V fjWv = which in the present case reduces to the single equation 

BB'{W^y + 2{BB" + B'^)Wj: + 2{BB" - B'^)Wg = 0. (2.13) 

We also calculated W!p explicitly and verified that indeed Wj^ = 0. 



3 String-Like Solutions 



Solving the equations becomes possible thanks to the following identity 



3 



-3/3/2 



,3/3/2 



/3" + 
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- -^'''e^f^/^ ( ^" + 



(3.1) 



which may be written as a total derivative in the following form: 



WI 



2e {€l3"f'^ < 



-3/3/2 



=3/3/2 

7W 



4 (e/3") 



//\3/2 -3/3/2 



(3.2) 



where e = ±1 according to whether /3" is positive or negative. The motivation to look for such 
identities comes from the fact that the components of Bach tensor may be expressed in terms of 
the Weyl tensor as in Eq. (1.4). Although it may seem simpler to get solutions in other gauges 
like the Mannheim gauge or in terms of the metric component H{r) defined above, we were able 
to integrate the equations in terms of B{r) (or /3(r)) only. 

Solving now the equation = is straightforward and quite simple. We will use both /3 
and B alternatively according to convenience. The first kind of solutions satisfies j3" = which 
is solved by 

B{r) = Boe^' (3.3) 

where k may be either positive or negative. 

The only other possibility is that the ratio [^\f7eP^'y jieT^)' in the second factor of (3.2) is 
a constant, say c, so we obtain after two integrations the following "mechanical" equation 

Cy^O ; \^{fi'f ^ec?e-^ = E , c = ^ (3.4) 

where a and E are integration constants. Note that the special case c = needs special care. 
A special family of solutions of this case is the above (3.3). Others will be considered later. 
Actually, it is easy to see from Eqs (2.11) and (3.2) that c = corresponds to conformally flat 
solutions with vanishing Weyl tensor. 

We therefore define an "effective potential" Veff{/3) which satisfies (/3')^/2 + Veff{/3) = E: 

, c = 0. 



i(/?0' + ^(a + ce- 



E 



VeffiP) 



£ 

3c^ 



(3.5) 



Fig. 1 presents the general behavior of Veff{f3) for e = +1. The curves for e = — 1 are 
just the same taken "upside-down". In terms of the metric function B the equation becomes 
{B')^/2 + Veff{B) = where now 



VeffiB) 



eg 
3c 



E)B^ + ea'^B + ^ + eac , c^O 



ea^B - EB^ 



c = 0. 



(3.6) 



The solutions may be expressed directly in terms of the metric function B, but part of the 
presentation is clearer in terms of p. So we will use both forms of solutions 

r{B) = J dB/^-2Veff{B) ; r(/3) = J d/3/^2 (£; - K//(/5)) (3.7) 



-1 1 2 3 4 5 6 



Figure 1: The effective potential V;//(/3) for e = 1, c = 1 and a = -1.25, -1, -0.75, -0.5, 0, 0.75. In 
order to identify tlie curves, note that the asymptotic value of 14// (/3) increases with a/c for e = 1. 

where the integration limits are determined by the boundary conditions. 

By inspection of the potential function one can conclude that there are 3 distinct types of 
solutions: one for e = +1 and two for e = — 1. It is easy to see that for e = +1 there are solutions 
only for values of E obeying E > /3c {or E > for c = 0). The effective potential 14//(/3) 
decreases monotonically with /3 and the solutions have a finite minimal value unless c = 0. On 
the other hand, for e = —1, E can have any real value. This case splits therefore according 
to whether E is above or below the maximal value of the effective potential Veff{/3) which is 
— a'^/3c (or 0). If £^ < —a^/3c, /3 or B is bounded from above. If E > —a^ /3c, the solutions can 
have the whole range < < oo. For future use we define the "energy excess" parameter C, by 



It is possible to express the solutions in terms of hyperelliptic integrals, but since the effective 
potential is monotonic, the possible behaviors of the solutions are limited so it is enough to 
introduce the following two kinds of real functions: 



3E/c^ = C^ + e{a/cf. 




u> I 



(3.8) 




< li < 1 



(3.9) 



and 




li > 



(3.10) 




max 



(3.11) 
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1 2 3 4 1 2 3 4 

(a) (b) 

Figure 2: Plots of the T-functions: (a) t[^-'(q;, m); the (+,— ) correspond to « > 1, u < 1 respectively, 
(b) T(±)(7,n). 

where Umax is the solution of (7 + =7^ + 1, that is: 

_ 1 

VI + 7^ -7 
with the cube root -y/x means — for j; < 0. 

Figure 2 shows the typical bchaAdor of the T-functions. Note that there are only three 
independent T-functions: T^ \a,u) and Tg ''(7,^^) are actually the same up to simple rescaling 
and translation: 

T^-\a,Q)-T^-\a,u) = ^T\-\^Xou) ; Co(«) = (1 + 3a + 3«2)V3 (3.13) 

In the special cases a = and 7 = (which correspond to a = 0) the T-functions get the 
elementary forms: 

t5+)(0,u) = ^ cosh-\u^/2) , t5")(0,u) = ^cos-^(n=^/2) 

T^+)(0,u) = -sinh-i(n3/2) ^ t(-^(0,«) = | sin-Htx'^/^). (3.14) 
3 3 

These correspond to the constant Ricci solutions that were mentioned in the previous section 
following Eq. (2.6). The fact that a = corresponds to constant Ricci solutions can be also 
inferred by noticing that in this case the potential function 14// (-B) of Eq. (3.6) reduces to that 
in Eq. (2.6). 

Two of the three kinds of the solutions of Eq. (3.4) can be expressed in terms of the Ti 
functions as 

y|^(.-..) = TP(>.|) (3.15) 



5 10 15 20 25 30 35 2 4 6 8 10 

(a) (b) 

Figure 3: Typical solutions for e — 1. Remember that L = B' . The other parameters are c = 1 and 
a — —1.25, —0.75, —0.5, 0, 0.25, 0.75. Note the non-monotonic dependence of the /3-slope on a. 



where the it corresponds to the e value. The e = +1 solutions are open with a minimum at 
B = Bq while the e = —1 ones are closed (that is bounded from above) for E < — a^/3c. 
Note that in this case B(r) can be extended by "reflection" with respect to the r = vq line, 
but L{r) = B' which vanishes there, reduces the domain of the solutions such that B(r) is 
monotonic. We will not elaborate on this kind of solutions further. 



n 1.0^ 




5 10 15 20 2 4 6 8 10 



(a) (b) 

Figure 4: (a) The Ricci components Rq (sohd) and (dashed), (b) The Weyl tensor component C. All 
plots for e = 1, c = 1 and a = —1.25, —0.75, —0.5, 0. These are 4 of the solutions shown in Fig. 3. The 
curves for the other two cases are similar, but have much more negative Ricci components. 
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Fig. 3 shows typical solutions for e = 1 with the boundary condition B{0) = 1 and -B'(O) = 
0. The components of the corresponding Ricci tensor are presented in Fig. 4. They are all 
compatible with the asymptotic condition (2.3) as should be the case. Indeed, it is easy to find 
from (2.5) that 

K = -Sc^C^ (3.16) 

where ( is the "energy excess" parameter defined above. Since > for the open solutions, all 
of them are asymptotically "negatively curved" (that is k < 0). 

A special member of this family is the AdS soliton which corresponds to a = and has 
constant Ricci components, as it solves Eq. (2.2) or (2.6). Note however, that unlike (A)dS 
space, it is not conformally flat, as is also obvious from the corresponding curve in Fig. 4b. 

The third type of solutions is with e = — 1 and E > —a^/3c, so they can have the whole 
range < B < oo. It may be written simply using the energy excess parameter C which is 
positive in this case too: 



.|C^/2(r 



TW(^,Ci?) 



(3.17) 



Fig. 5 shows typical solutions of this kind where we choose -B(O) = 0. Fig. 6 has the 
components of the corresponding Ricci tensor which shows a singularity on the axis. However, 
these solutions may still be physically relevant as exterior solutions of appropriate cylindrically- 
symmetric sources. Wc took the same values as in Fig. 3 for the parameters a and c as well 
as for the "energy excess" i^. We also chose to present sinh^^(i?) instead of B in order to 
cover both the small B and large B regions. The exponential increase of B(r) away from the 
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Figure 5: Typical open solutions for e = — 1. Remember that L = B' . The other parameters are c = 1 
and a = —1.25, —0.75, —0.5, 0, 0.25, 0.75. Note again the non-monotonic dependence of the slopes on o. 
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0123456" 1 2 3 4 5 

(a) (b) 

Figure 6: (a) The Ricci components Rq (solid) and (dastied). (b) Tlie Weyl tensor component C. All 
plots for e = —1, c — 1 and a = —1.25, —0.75, —0.5, 0, 0.25. These are 5 of the solutions shown in Fig. 
5. The a — 0.75 curves are similar to the a = 0.25 ones, but have much more negative Ricci components. 

axis is obvious. An exception from the generic exponential behavior is the minimal £'-value for 
unbounded solutions which correspond to C = 0. In this case the asymptotic behavior is i? ~ r^. 
The asymptotic behavior of the Ricci components is given again by (2.3) with (3.16) so in both 
cases K < 0. The C = solution is asymptotically Ricci flat. 

Note that again the special case a = allows for explicit (constant Ricci) solutions in terms 
of elementary functions - see (3.14). For the "last unbounded trajectory", which corresponds to 
C = (with e = —1) we find B = [^y3/2\c\{r — vq)]'^^^. This solution is not only asymptotically 
Ricci flat, but has Rfj_u = identically and is actually the well-known "Kasner" string-like 
solution of ordinary GR. Another special case which requires anyhow a separate treatment is 
c = 0. In this case we can just integrate Eq. (3.7) directly to get the general solutions 

B = B„c„sh^(M^) , , = 1 (3.18) 
B = B„cos=(!ii^) . . = -1, E<0 

instead of (3.15), while instead of (3.17) we have 

i2 



B = — sinh^ (r - ro) j , e = -1 , E>0 (3.19) 

The special case E = can be obtained from the limit of this equation to give a behavior of 



B = a'^{r -rof/2 , e = -1 , E = 
11 



(3.20) 



This family is supplemented by the exponential solution (3.3) encountered already at the begin- 
ning. 

All these c = solutions have non-trivial curvatures, but as mentioned above they are 
conformally flat. However, they are not trivial even within the framework of CG, since they are 
either conformal to Minkowski space-time with a conic angular deficit like (3.20), or to a certain 
region of it like the others. The solution (3.20) is therefore the analog in this gauge of CG to 
the conical string-like solution of GR. 

Finally we note a possible generalization of this work to the scalar-tensor extension of CG 
which allows for a breakdown of the conformal symmetry. Some cylindrically-symmetric solu- 
tions of this kind were obtained numerically in a recent study [12], but considering the results 
presented in this section, it may turn out possible to obtain explicit exact solutions in the gauge 
where the scalar field is constant. The metric tensor will have now two independent compo- 
nents, but the second may be taken as a conformal factor (as in ref [7]) which will simplify the 
equations so that analytical treatment will be possible. 



4 Null Geodesies 

We proceed here to examine further the nature of these solutions by studying their geodesies. 
The geodesic equations in the metric diag{B^, —1, — L^, —B^) are easily integrated to give 

z = k , L'^ip/B^ =i , f^ + {k'^- 1)B^ + fB'^/L'^ = -nB^/S'^ (4.1) 

where the coordinates are functions of time t and k < 1, i and £ are constants of the motion. 
The parameter fi = 0,1 distinguishes between lightlike or timclikc geodesies respectively. Since 
only null geodesies have an invariant meaning, we will concentrate on them. 

One can analyze numerically the trajectories by the effective potential for the r-motion which 
satisfies r^/2 + U,.ff{r) = 0. However, analytic treatment is possible if we "compromise" on a 
parametric representation of the solutions r[t) by {r{B),t{B)). We will not expect any trouble 
with that since we concentrate in the open solutions where B{r) is monotonic. The r-equation 
wiU be replaced by B^/2 + Weff{B) = with 



WeffiB) 



B^ + t^{\-k^)[o?B^ ^aB+'^\B (4.2) 



where a = ajc. This is the potential for the general case with c 7^ 0. The c = case needs 
as usual a special treatment and is given by (4.8) below. Actually, since there exist explicit 
solutions in this case one may obtain also Ueff{r) explicitly. It is straightforward to analyze 
the different trajectory types from Weff{B) and especially its zeroes. The solutions may be 
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bounded (in B or in r) where Weff{B) < in a finite interval or unbounded (or open) if this 
interval extends to infinity. The process is simplified by noticing that W(,ff{B) has only one 
zero in addition to the possible one at B = 0. 

For e = 1 all solutions are open, namely Bmin < B < oo where Bmin depends on the 
constants of the motion k and I and on the other geometrical parameters a, c and E (see 
below). A necessary condition for solutions to exist (for e = 1) is that the coefficient of the B^ 
term in Weff{B) which is proportional (with the positive coefficient c^(l — fc^)) to 

will be negative, otherwise there will be no zeroes of Weff{B) for -B > 1 which is the domain of 
B in this case. This imposes an upper bound on the combination ^^/(l — fe^): 

<2[E- —— (4.4) 



1- P 



For e = — 1 there are two cases: open solutions which now extend for < < oo occur if Eq. 
(4.4) is satisfied. In the complementary case we have bounded solutions for < i? < B^ax where 
Bmax depends on the characteristic parameters. Bmax and also Bmin can be written together as 

2 

Bminlmax = ^ ("^ + ^^^/a^ - 1 - ^(e^/a^ _ 1)2^ (4.5) 
with the corresponding value of e. 




(a) (b) 

Figure 7: Effective potential (divided by (1 — fc^)/2) curves for null geodesies in the two types of 
metrics with c = 1 and a = -0.5. (a) e = 1, ^ = -0.15, -0.13, -0.11, -0.10, -0.09 ; (b) e = -1, 
^ = —0.2, 0, 0.07, 0.1, 0.15. In both cases the higher curves correspond to the larger values of ^. 
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Plots of Ueff{r) which present ah these properties are shown in Fig 7. For e = 1 only ^ < 
curves are shown since otherwise U(.ff{r) is always positive and no solutions exist. Note that 
each curve in this figure corresponds to a single null geodesic (up to ^-translation) defined by k 
and I since only vanishing "effective energy" is allowed. 

In order to find out about the shape of the geodesies, one may solve for r{(p) or equivalently 
for the parametric representation {r{B),ip{B)). It turns out to be simpler to transform to 
V = 1/B where we find the following first order equation for light-like trajectories, or more 
accurately, their projection on the (r, ip) plane: 



1 fdvV 



1 



2(1 -k^ 



E- 



ec 



(a + vf 



E- 



ec 



[a + v) 







(4.6) 



£2 \~ 3 

The projections of some typical trajectories appear in Figs. 8-9. They demonstrate explicitly 
the features obtained from the general discussion above. The e = 1 solutions are all open with a 
"periastron" chosen to be always at ip = 0. Several windings are possible according to the values 
of the parameters. The e = — 1 solutions may be either open with < r < cxd, or bounded with 
< r < 

^max- Again, several windings are possible. 
Finally we discuss shortly some special cases. The first is a = 0. Generally null geodesies in 
this case are just slightly different from the adjacent solutions except when also ^ = for e = — 1. 
This is the only case when Ueff{r) vanishes asymptotically to allow open geodesies with r{t) — )• 
as t — )• CXD. The resulting orbits are spirals which become denser as r increases. Actually, the 
effective potential U(,ff{r) can be written explicitly in terms of elementary functions, and the 



\ e=+l';c=l y 
_ \ a=-0.5 \ 






^^^^^^^^^ 


(=-0.01 ^\ 
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-20 -10 10 



(b) 



Figure 8: Projection in the {r,(p) plane of null trajectories with e = +1, c = 1, a = —0.5. (a) Solutions 
with (^-values from ^ = —0.07 to f = —0.01 in even steps; (b) Three more winding solutions with 
C = -0.11, -0.10, -0.09. 
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(a) (b) 

Figure 9: Projection in the {r,(p) plane of null trajectories with e = — 1 , c = 1, a = —0.5. (a) Open 
solutions with ^-values from ^ = —0.2 to ^ = in even steps; (b) Bounded solutions with ^-values from 
^ = 0.2 to ^ = 1 in even steps. 

typical null geodesies with r(0) = and tp{Q) = are 

r{t) = Y^^sinh-i(?) ^ ^ IL^jZIZ i- ^ (tan-\2i-V3)+tan-\2t + x/S)) 

-- tan-^ i + V31og -y=^ (4.7) 

3 ^yp + ^t + lj\^ ^ 

where t = ^y{l-k^)/6 ct. 

The case c = also allows to obtain explicit expressions for Ueff{r), but it is easier to solve 
for B{t) using the effective potential of Eq. (4.2) which now simplifies to 

<7^(i?) = - (1 - k')E^ + 6(1 - k')a'B^ (4.8) 

The transformation v = 1/B gives a mechanical equation with a linear potential which gives 
easily the following solutions in the various cases. Two of them can be written as 

B- 2 — + 2a^J 

where we chose t = to be the time which corresponds to the minimal or maximal radial 
distances for e=lore = — 1 respectively. For e = 1, t should be further restricted such 
that 1/B is non-negative. Note that the relative sizes of the two quantities in the bracket of 
the second term are such that it is always positive. The third kind of geodesies which have 
< i? < oo exist only for e = — 1 and are given by 

V +V'^{l-k^)E-£H (4.10) 
B 2 
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In order to get the explicit r{t) dependence of these geodesies one needs the corresponding B{r) 
solutions for this case - Eqs. (3.19)-(3.19). We will not present here this last step. The explicit 
solution ip{t) which we skip too may be obtained by direct integration of (see (4.1)) 

^ = ^^ = (411) 

dt L2 2{E-ea?/B) ^ > 

with l/B{t) obtained from (4.9)-(4.10) above. 



5 Conclusion 

We have analyzed the vacuum field equations of CG in the static cylindrically-symmetric case and 
found all solutions explicitly. In some cases the solutions are expressed in terms of elementary 
functions. 

There are three kinds of solutions: two open ones and one closed. The open ones split into a 
regular family (with a possible conic singularity on the axis) which contains the AdS soliton and 
a family of singular solutions where 500 vanishes on the axis and the Ricci and Weyl components 
diverge there. These two families of open solutions have generically asymptotically vanishing 
Weyl tensor. Some special cases are conformal to conical space-times. 

We have further analyzed the null geodesies in the two families of the open spacetimes. 
We found that the singular spacetimes support two kinds of null geodesies: open ones (that is 
< r < 00) which are usually spirals when the angular momentum is below a certain maximal 
value, and bounded ones when the angular momentum is above the critical value. In the regular 
spacetimes only open geodesies exist for the same range of angular momentum. Unless the 
angular momentum vanishes, the orbits in this case avoid r = 0. Above the critical value no 
geodesies exist at all. 

A possible future application of these results is to use these families of exact solutions to 
analyze light bending in the vicinity of localized linear sources in CG. The analogous problem 
of light bending in spherically-symmetric gravitational field was only recently settled [19]. 
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